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A PROCEDURE  FOR  EVALUATING  FRACTURE  DYNAMIC  PARAMETERS  FROM  CRACK 
VELOCITY  MEASUREMENTS 


A.  S.  Kobayas’  and  S.  Mall 

Department  of  Mechanical  Engineering,  University  of  Washington 
ABSTRACT 

A combined  numerical  and  experimental  procedure  for  evaluating  some  of 
the  fracture  dynamic  parameters  which  govern  the  crack  run-arrest  response 
in  a fracturing  plate  are  discussed.  A dynamic  finite  element  code  is  used 
to  compute  the  dynamic  stress  intensity  factor  and  dynamic  energy  release 
rate  associated  with  a propagating  crack  which  is  driven  by  the  experimen- 
tally determined  crack  velocity.  Numerical  results  generated  by  the  de- 
veloped procedure  are  then  compared  with  dynamic  stress  intensity  factors 
obtained  through  dynamic  photoelastic  analyses  of  fracturing  Homalite-100 
plates.  Two  edge-cracked  specimens  with  fixed  edge  displacement  loadings 
and  two  wedge-loaded  double  cantilever  beam  specimens  were  considered  in 
this  comparative  study.  Good  agreements  were  obtained  between  the  results 
obtained  by  the  developed  numerical-experimental  procedure  and  dynamic 
photoelasticity. 

INTRODUCTION 

The  three  approaches  currently  in  use  in  fracture  dynamic  studies  are: 
to  relate  the  experimentally  determined  crack  velocities  with  static  frac- 
ture parameters  [1];  to  relate  experimentally  determined  crack  velocities 
with  those  obtained  from  an  analytical  dynamic  model  withopt  or  with  postu- 
lated dynamic  fracture  toughness  [2,3,4];  and  to  determine  experimentally 
Che  dynamic  state  of  stress  in  fracturing  polymeric  materials  [5,6,7],  In 
this  paper,  a fourth  procedure,  which  utilizes  the  versatility  of  dynamic 
finite  element  method  (FEM)  to  extract  fracture  dynamic  parameters  from  the 
most  commonly  measured  quantity  of  crack  velocity  in  practical  structural 
material  is  described.  The  procedure  is  then  used  to  indicate  the  errors 
involved  in  using  static  analysis  to  Interpret  dynamic  results. 

DYNAMIC  FINITE  ELEMENT  CODE 

The  dynamic  finite  element  method  (FEM)  which  has  evolved  since  its  ini- 
tial use  [8]  in  this  combined  numerical-experimental  procedure  is  the  FEM 
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code  HONDO  [9)  with  artificial  viscosity  to  reduce  the  keystoning  effect  in 
the  finite  elements  along  the  opening  crack  surfaces.  In  this  dynamic  FEM 
procedure,  crack  extension  is  modeled  by  discontinuous  Jumps  of  the  crack 
tip  from  one  finite  element  node  to  its  adjacent  node  and  the  time-averaged 
. rack  tip  displacement  is  equated  to  the  average  crack  tip  velocity  at  the 
adjacent  node.  The  average  work  necessary  to  open  the  newly  created  incre- 
me.nt  ol  crack  surface  is  then  used  to  compute  the  surface  energy  dissipation 
rate  which  is  equated  to  the  dynamic  energy  release  rate  at  the  crack  tip 
node  just  prior  to  the  subsequent  discrete  movement  of  the  crack  tip  [8]. 

An  .iccuracy  check  of  this  direct  procedure  for  computing  the  dynamic  energy 
release  rate  was  established  by  a comparative  study  with  Baker's  solution 
IIG]  where  the  numerical  and  theoretical  results  agreed  within  1 percent  of 
each  ocher  [8j. 

More  recently  another  check  on  the  dynamic  finite  element  algorithm  was 
made  by  comparing  the  crack  opening  displacements  (COD)  of  a constant  velo- 
city crack  against  the  theoretical  solution  of  Broberg  [11].  In  tliis  numer- 
ical study  involving  a fracturing  steel  plate,  an  artificial  keystone  vis- 
cosity of  B3=0.2  was  used  to  damp  out  the  keystoning  effect  prevalent  in 
previous  analyses.  The  numerically  determined  time-averaged  COD  at  every 
other  node  away  from  the  moving  crack  tip  was  found  to  be  in  excellent  agree- 
ment with  Broberg' s result  for  a crack  velocity  of  c/cp  = 0.076  where  c and 
c are  the  crack  and  dilatational  wave  velocities,  respectively  [12].  It 
was  also  shown  that  at  this  low  crack  velocity,  the  stress  intensity  factor 
computed  by  static  near  field  solution  was  only  2.2  percent  higher  than  that 
computed  by  the  dynamic  near  field  solution  when  dynamic  finite  element  COD 
values  were  used  to  compute  the  dynamic  stress  intensity  factor. 

In  order  to  examine  further  the  effectiveness  of  the  above  COD  technique 
in  dynamic  finite  element  analysis  of  fracturing  Homalite-100  plates,  the 
weJl  analyzed  fracturing  dynamic  photoelastic  specimens  B2  and  B13  [8]  were 
reanalyzed  with  various  artificial  viscosities.  The  finite  element  break- 
down and  the  crack  velocities  used  in  this  new  study  are  identical  to  those 
used  in  Reference  [8].  Figure  1 shows  the  crack  opening  displacements  at 
three  crack  lengths  in  Specimen  B2  with  three  artificial  viscosities.  The 
simill  artificial  viscosity  of  B3  = 0.01  and  0.1  were  not  very  effective  in 
removing  keystoning  but  the  prominent  fluctuations  in  COD  were  substantially 
suppressed  with  B3  = 0.5.  It  should  also  be  noted  that  a significant  change 
was  made  in  the  numerical  algorithm  where  the  initial  residual  surface  trac- 
tions along  crack  propagation  were  computed  by  using  a static  modulus  of 
elasticity  of  Es  = 540  ksi  while  all  stress  wave  propagation  induced  by  the 
running  crack  were  computed  with  a dynamic  modulus  of  elasticity  of  Ep  = 

675  ksi.  This  combined  use  of  static  and  dynamic  moduli  of  eiastic^'y  in 
dynamic  finite  element  analysis,  initially  proposed  by  Gehlen  [1j],  an 
attempt  to  model  the  strain  rate  sensitivity  of  the  modulus  of  elasticity 
of  Homalite-100  material.  Such  strain  rate  sensitivity  did  not  exist  in  the 
fracture  dynamic  analysis  of  steel  tapered  DCB  specimen  [12]  which  is  easiei 
to  analyze  by  dynamic  finite  element  analysis  than  the  photoelaslic  speci- 
mens. Our  continuing  Interest  in  using  dynamic  photoelastlcity  to  study 
fracture  dynamics  despite  the  added  complexity  of  strain  rate  sensitivity 


Fig.  (1)  - Effect  of  artificial  keystone  viscosity  on  COD  in 
Test  No.  B2 . 


stems  from  the  fact  that  no  other  experimental  technique  can  provide  an  ac- 
. I curate  near  field  state  of  stress  in  the  vicinity  of  a running  crack.  Dy- 

namic photoelasticity  is  one  of  few  optical  techniques  which  can  provide 
experimentally  determined  dynamic  stress  intensity  factor  with  which  our 
numerical  results  could  be  compared.  As  will  be  shown  later,  such  dynamic 
stress  intensity  factors  computed  through  the  combined  use  of  static  and 
dynamic  moduli  of  elasticity  agreed  well  with  the  experimental  results  under 
the  state  of  plane  strain.  The  latter  state  of  plane  strain  for  thicker 
Homalite-100  plates  is  considered  to  be  a better  modeling  of  the  cleavage 
fractured  surface  observed  in  the  fractured  specimens. 

Figure  2 shows  the  variations  in  dynamic  stress  intensity  factors  ob- 
tained numerically  by  the  energy  release  rate  method  [8]  and  COD  method  [12] 
as  well  as  their  experimental  counterpart  determined  by  dynamic  photoelas- 
ticity for  Specimen  B2.  While  the  two  numerical  algorithms  of  computing 
dynamic  stress  intensity  factors  do  not  yield  substantially  different  re- 
sults, the  more  pronounced  fluctuation  in  dynamic  stress  intensity  factors 
computed  by  the  COD  method  is  noted.  Thus  this  comparison  favors  the  ener- 
gy release  rate  procedure. 

Likewise  comparison  is  shown  in  Figure  3 for  Specimen  B13  in  which  the 
crack  arrested.  While  the  prominent  peak  in  the  experimentally  determined 
dynamic  stress  Intensity  factor  was  not  observed  in  the  two  numerical  values, 
the  experimental  and  numerical  results  are  otherwise  in  good  agreement  with 
each  other.  Again  the  numerical  results  obtained  by  energy  release  rate 
procedure  appears  to  be  in  slightly  better  agreement  with  the  experimental 
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Fig.  (2)  -Dynamic  stress  Inten- 
sity factors  in  Test 
No.  B2,  plane  strain 
analysis 


Fig.  (3)  - Dynamic  stress  inten- 
sity factors  in  Test 
No.  BIJ,  plane  strain 
analysis 


As  a result  of  these  two  comparative  studies  using  the  well  analyzed 
specimens  B2  and  B13,  it  was  concluded  that  the  energy  release  rate  proce- 
dure provided  slightly  better  numerical  result.s  of  the  fracturing  Homalile- 
100  plates.  This  conclusion  unfortunately  differs  with  that  in  Reference 
[12]  and  is  perhaps  indicative  of  the  larger  keystoning  effect  in  the  Homa- 
llte-100  plates  due  to  its  much  smaller  modulus  of  elasticity. 

WEDGE-LOADED  DCB  SPECIMEN 

The  above  dynamic  finite  element  algorithm  was  then  used  to  analyze  the 
dynamic  crack  arrest  results  of  two  wedge-loaded  double  cantilever  beam 
(DCB)  specimens  machined  from  Homallte-100  plates  of  1/2  inch  nominal  thick- 
ness [14].  The  dynamic  stress  intensity  factors  of  these  two  fracturing 
tapered  DCB  specimens  were  determined  by  a different  data  reduction  scheme 
of  dynamic  photoelasticity  results  than  that  used  by  the  authors  (7).  The 
static  and  dynamic  material  properties  as  well  as  the  crack  position  versus 
time  relations  reported  in  Reference  [14]  were  used  as  input  conditions  to 
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our  dynamic  finite  element  analysis. 

Figure  4 shows  the  finite  element  breakdown  of  one  of  the  two  tapered 
DCB  specimens.  Figures  5(a)  and  (b)  show  the  crack  position  versus  time 
relations  which  were  used  to  propagate  the  crack  tip  intermittently  along 


total  NODE:  390 


Fig-  (^)  “ Finite  element  breakdown  of  wedge-loaded  C-DCB 
model  No.  7 [14] 


Fig.  5(a) 


Crack  tip  position  versus  time  in  C-DCB  model 
No.  6 [14] 
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Fig.  5(b)  - Crack  tip  position  versus  time  in  C-DCB  model 
No.  7 [14] 


tbe  finite  element  nodes  at  prescribed  time  intervals.  Fixed  pin  displace- 
ments were  prescribed  during  the  entire  fracture  process  which  were  459  and 
100  microseconds  for  Specimens  C-DCB  Model  Nos.  6 and  7,  respectively.  This 
assumption  of  fixed  grip  condition  has  been  the  source  of  discussion  since 
one  of  the  authors  presented  a dynamic  finite  element  analysis  of  the  crack 
propagation  and  arrest  in  a tapered  DCB  specimen  using  measured  crack  velo- 
cities [15].  The  static  stress  Intensity  factor  at  crack  arrest  computed 
under  tile  assumption  of  fixed  grip  condition  was  approximately  40  percent 
lower  than  tlie  corresponding  value  computed  under  the  assumption  of  varia- 
ble load  condition  of  Reference  [16],  thus  Indicating  the  seasitivity  of 
tile  arrest  stress  intensity  factor  to  the  prescribed  boundary  condition 
during  crack  propagation  in  such  small  crack  arrest  specimen.  Subsequent 
analysis  of  another  tapered  DCB  specimen  under  the  two  different  loading 
conditions  of  our  fixed  grip  loading  and  the  variable  loading  condition  pre- 
scribed in  References  [16]  and  [17]  showed  that  while  the  dynamic  strains 
computed  under  the  former  condition  agreed  well  with  the  three  strain  gage 
measurements  [lo],  the  corresponding  dynamic  strains  computed  under  the 
assumption  of  variable  loading  differed  considerably  with  experimental  re- 
sults [12].  As  a result  of  this  analysis  [12],  the  fixed  grip  condition  is 
believed  to  be  a valid  assumption  in  the  series  of  experiments  reported  in 
Reference  [16].  This  conclusion  also  Indicates  that  the  crack  arrest  stress 
intensity  factor  determined  by  the  procedure  described  in  References  [16] 
and  [17]  could  grossly  overestimate  the  crack  arrest  potential  of  the  mater- 
ial tested . 

Unfortunately,  the  excellent  agreements  between  the  crack  arrest  stress 
intensity  factors  determined  by  simulations  of  the  crack  arrest  experiments 
of  References  [16]  and  [17]  using  Horaali te-100  specimens  [14]  and  the  crack 
arrest  stress  intensity  factors  determined  independently  by  dynamic  photo- 
elasticity [ 5]  have  been  used  as  experimental  evidence  for  justifying  the 
variable  pin  loading  in  References  [16]  and  (17).  A cursory  study,  however. 


shows  considerable  differences  between  the  relative  compliances  between 
the  loading  fixture  and  the  Homalite-100  specimen  in  the  dynamic  photoelas- 
ticity experiments  and  those  of  References  [16]  and  [17].  As  a result,  one 
can  conclude  qualitatively  that  the  variable  pin  loads  measured  in  the  ri- 
gid loading  fixture  of  the  simulated  experiments  using  Homalite-100  speci- 
mens should  be  closer  to  the  pin  loads  obtained  under  fixed  grip  condition 
while  such  condition  cannot  be  realized  *n  the  actual  experiments  using 
steel  specimen.  An  experimental  check  in  these  relative  compliances  can  be 
easily  made  to  verify  such  hypothesis.  Another  procedure  is  to  analyze 
numerically  the  fracture  dynamic  response  of  the  Homalite-100  specimens 
following  the  procedure  described  in  Reference  [15]  and  then  identify  the 
differences  or  similarities  between  these  results  with  those  of  Reference 
[12],  Dynamic  finite  element  analyses  under  the  fixed  grip  condition  should 
thus  provide  this  insight  into  the  controversy  surrounding  the  exact  boun- 
dary conditions  on  the  tapered  DCB  specimens  used  in  crack  arrest  experi- 
ments . 

Figure  6 shows  the  dynamic  and  static  stress  intensity  factors  obtained 
by  static  and  dynamic  finite  element  analyses  under  the  assumption  of  fixed 
grip  condition  for  a Homalite-100  C-DCB  Model  No.  6 [14].  Also  shown  are 


CRACK  EXTENSION,  0-a„,  in 


Fig.  (6)  - Stress  intensity  factors  in  a wedge-loaded  contoured 
DCB  specimen  (Model  No.  6).  Load  at  fracture  initia- 
tion = 151  lb  [14]. 


Ltu  dyrainic  stress  intensity  tactors  of  tills  specimen  obtained  by  dynamic 
piiotoeiast  ic  ity  ll^J.  While  tlie  nnnierically  obtained  dynamic  stress  inten- 
sity factors  sliow  some  sm.ill  oscillations,  Ctie  two  results  are  j’enerally 
in  qualitative  ap,reenient  witli  eacli  otlier  and  in  particular  are  in  excellent 
agreement  at  crack  arrest. 


the  crack  opening  displacements  (COD) 
, tile  keystoning  effect  continues  to 


Crack  opening  displacement  (COD)  at  different  times 
in  C-DCB  Model  No.  6 [14] 


increase  in  crack  lengtli  despite  the  high  artificial  viscosity  of  B3  = 0,5 
The  pronounced  smaller  dynamic  crack  opening  displacements,  wliich  was  also 
observed  in  tlie  metallic  specimens  [12,15],  in  this  cantilever  beam  type 
specimen  are  indications  of  the  delayed  response  of  the  propagating  crack 
tip  to  the  applied  load  at  the  loading  pin.  Such  delay  response  further 
verifies  our  assumed  fixed  grip  loading  condition  during  crack  propagation 
and  arrest. 


Figure  8 sliows  the  dynamic  and  static  stress  intensity  factors  in  Speci- 
men C-DCB  Model  Mo.  7.  Again,  relatively  good  agreement,  particularly  at 
crack  arrest,  between  ttie  experimental  and  numerical  results  is  noted.  The 
rapid  small  oscillations  in  dynamic  stress  intensity  factors,  which  were 
noted  in  Figure  6,  are  absent  in  Figure  8.  Our  past  e.xperiences  in  numeri- 
cal fracture  dynamic  analyses  [A,  12]  indicate  that  smootlied  crack  veloci- 
ties generally  result  in  oscillations  in  dynamic  stress  intensity  factors 
and  vice  versa.  Since  the  test  results  of  Model  C-DCB  No.  6 recorded  9 
crack  positions  for  a crack  extension  of  J.41  in.  while  7 crack  positions  are 
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Fig.  (8)  - Stress  intensity  factors  in  a wedge-loaded  contoured 
DCB  specimen  (Model  No.  7).  Load  at  fracture  initia- 
tion = 89  lb  [14]. 


recorded  for  a crack  extension  of  0.90  in.  for  Model  7,  one  would  expect  the 
crack  position  versus  time  relation  for  the  test  of  Model  6 to  represent  a 
smoother  time  average  position  instead  of  a more  precise  crack  position  ver- 
sus time  relation  necessary  to  generate  a smoothly  varying  dynamic  stress 
intensity  as  was  the  case  of  Model  No.  7 test.  Thus  Model  No.  7 should 
yield  smoother  variations  in  numerically  computed  dynamic  stress  intensity 
factor  than  Model  No.  6. 

The  above  good  agreement  between  the  variations  in  dynamic  stress  inten- 
sity factors  with  crack  propagation  obtained  numerically  and  experimentally 
verifies  the  validity  of  the  fixed  grip  condition  under  which  the  dynamic 
stress  intensity  factors  were  computed.  Although  the  numerical  results  pro- 
vide the  variations  in  pin  loads  with  respect  to  crack  propagation  as  shown 
in  Figure  9,  corresponding  experimental  results  were  not  available  for 
direct  comparison.  Thus  experimentally  measured  pin  load  from  a separate 
test  [14]  is  shown  in  Figure  9 for  qualitative  comparison.  The  two-fold 
differences  between  ringing  frequencies  of  the  numerical  and  experimental 
pin  loads  could  be  due  to  unavoidable  compliance  of  the  loading  frame  and 


T iVE  . ^sec 

Fig.  (9)  - Pin  load  versus  time  in  wedge-loaded  C-DCB  specimens 

load  cell  system  in  the  actual  experimental  setup. 

DISCUSSION 


f An  internal  accuracy  check  of  the  dynamic  finite  element  algorithm  was 

: made  by  computing  the  instantaneous  energy  balance  of  the  entire  system, 

' Typical  results  for  four  crack  lengths  are  shown  in  Tables  1 and  2.  It  is 

; interesting  to  note  that  the  accuracy  of  this  energy  balance  is  somewhat 

I lower  for  Model  No.  6 specimen  which  also  snowed  larger  oscillation  in 

I . dynamic  stress  intensity  factors  with  crack  propagation. 

I ’ I CONCLUSIONS 


A dynamic  finite  element  algorithm  has  been  developed  for  computing  the 
dynamic  stress  intensity  factors  from  experimentally  determined  crack  posi- 
tion versus  the  relation  of  a propagating  and  arresting  crack.  Accuracy  of 
the  developed  numerical  procedure  was  checked  by  comparing  the  numerically 
determined  stress  intensity  factors  with  those  determined  in  foui  dynamic 
plio toelas tic  experiments. 

The  numerically  and  experimentally  determined  dynamic  stress  intensity 
factors  of  two  wedge-loaded  DCB  specimens  indicate  that  the  crack  propagates 
under  a fixed  grip  condition. 
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TABLE  1 - ENERGY  BALANCE  IN  C-DCB  MODEL  NO.  6 


Crack 

Extension 

W 

u 

-I 

K 

F 

(U  + K + F)/W  - 1 

1.0 

3.201 

2.221 

0.223 

0.534 

- 0.070 

2.0 

3.201 

1.956 

0.074 

0.879 

- 0.091  ; 

3.0 

3.201 

1.748 

0.038 

1.114 

- 0.094 

3.5* 

3.201 

1.629 

0.055 

1.154 

- 0.113 

J 

TABLE  2 - ENERGY  BALANCE  IN  C-DCB  MODEL 

NO.  7 

Crack 

Extension 

W 

U 

K 

F 

(U  + K + F)/W  - 1 

0.25 

1.084 

0.079 

- 0.033 

0.50 

1.084 

0.867 

0.034 

0.116 

- 0.067 

0.75 

1.084 

0.855 

0.015 

0.155 

- 0.059 

0.875* 

1.084 

0.836 

0.005 

0.181 

- 0.062 

* Crack  Arrest 


W:  External  Work  Unit:  Pound-in. 

U:  Strain  Energy 

K:  Kinetic  Energy  fC  i 

F;  Fracture  Energy  j dc 

^ o 

c:  Length  of  Crack  Extension 
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A combined  numerical  and  experimental  procedure  for  evaluating  some  of  the  I 

fracture  dynamic  parameters  which  govern  tTie' crack  run-'arrest  response  in  a | 

fracturing  plate  are  discussed.  A dynamic  finite  element  code  Is  used  to  com- 1 
pute  the  dynamic  stress  Intensity  factor  and  dynamic  energy  release  rate  ! 

associated  with  a propagating  crack  which  Is  driven  by  the  experimentally  I 

determined  crack  velocity.  Numerical  results  generated  by  the  developed  pro- 
cedure are  then  compared  with  dynamic  stress  intensity  factors  obtained  ttiroug: 
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ABSTRACT  (Continued) 


edge-cracked  specimens  with  fixed  edge  displacement  loadings  and  two  wedge-loaded 
double  cantilever  beam  specimens  were  considered  in  this  comparative  study. 

Cood  agreements  were  obtained  between  the  results  obtained  by  the  developed 
numerical-experimental  procedure  and  dynamic  photoelasticity. 


